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Abstract 



We consider the Gibbons- Hawking metric for a three-dimensional periodic array of multi- 
r> ' Taub-NUT centers, containing not only centers with a positive NUT charge but also ones 

5t , with a negative NUT charge. The latter are regarded as representing the asymptotic form 

of the Atiyah-Hitchin metric. The periodic arrays of Taub-NUT centers have close parallels 
with ionic crystals, where the Gibbons-Hawking potential plays the role of the Coulomb 
static potential of the ions, and are similarly classified according to their space groups. After 
a periodic identification and a Z2 projection, the array is transformed by T-duality to a 
system of NS5-branes with the SU(2) structure, and a further standard embedding yields, 
though singular, a half-BPS heterotic 5-brane background with warped compact transverse 
dimensions. A discussion is given of the possibility of probing the singular geometry by 
two-dimensional gauge theories. 
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1 Introduction 



Branes and singularities are essential ingredients in modern string theory. On one hand, 
branes are non-perturbative objects whose existence is required for consistency of various 
dualities. On the other hand, singularities are sources from which the effects of such objects 
emerge by allowing them to wrap around collapsing cycles. Branes and/or singularities are 
also required in order to realize inflation in string theory; under some mild assumptions, 
it was shown that smooth compactifications cannot allow an accelerating solution [1]. One 
way to circumvent this is to consider warped compactification with branes or singularities. 
There is a growing literature on attempts to realize a Randall- Sundrum- like framework [2] 
in string theory by using D-branes [21 HI E] . 

It has been known that, in order to achieve warped compactification in supergravity, some 
of the constituent branes need to have negative tension [3] . (This fact was shown in type IIB 
theory but can be generalized to other string theories.) In type II theories, an orientifold has 
this property. In heterotic string theory, however, one cannot define orientifolds as one does 
in type II theories, since, while one needs to gauge world-sheet parity to define orientifolds, 
heterotic theory is left-right asymmetric and hence does not possess such a symmetry to be 
gauged. Although even in the absence of the worldsheet definition one can define orientifolds 
by simply assigning Z2 actions on various fields [B], it is unclear what kind of geometry is 
produced there. The back- reaction problem of orientifolds was considered in [7]. 

In this paper, we propose a novel, interesting set up for warped compactification of heterotic 
(as well as type II) string theory to six dimensions, based on a simple idea: We consider 
the Gibbons-Hawking metric [8] for a three-dimensional periodic array of multi-Taub-NUT 
centers consisting of not only centers with a positive NUT charge but also ones having a 
negative NUT charge. This negative- charge metric is regarded as representing the asymptotic 
form of the Atiyah-Hitchin metric; as is well known, the Atiyah-Hitchin metric [S] is a smooth 
self-dual metric on a Euclidean four-manifold, and it looks like a negative-charge Taub-NUT 
for an observer who is far enough from the bolt locus. Such a periodic array of infinite Taub- 
NUT centers was already discussed in [TOl |Tl]. The idea of introducing an Atiyah-Hitchin 
space in the Gibbons-Hawking metric was used in [12] to understand how an S0{2n) gauge 
symmetry enhancement occurs due to D-branes. By T-duality, the periodic array of Taub- 
NUT centers is converted to a periodic array of NS5-branes. An ordinary (positive-charge) 
Taub-NUT is T-dualized to a (smeared) NS5-brane, and the negative-charge Taub-NUT is 
mapped to a (smeared) negative-tension brane. The resulting T-dual system has the SU{2) 
structure (hyper-Kahler with torsion), and preserves one half the supersymmetries. We then 
embed the generalized spin connection into gauge connection to obtain a 5-brane system 



with a compact four-dimensional transverse space in heterotic string theory. 

As is known, although the negative-NUT approximation of the Atiyah-Hitchin metric is 
good at large distances compared to the NUT-charge scale, near the center the approximate 
metric becomes singular on a two-dimensional sphere (in the three-dimensional space where 
the periodic array is defined) of radius equal to the NUT charge. Therefore, on the T-dual 
side, there are two kinds of regions where the supergravity analysis is not reliable: One is 
the region near a smeared NS5-brane, the familiar one where the geometry is an infinite 
throat and the dilaton blows up as one approaches the core. The other is the region near the 
brane dual to the nagative-charge Taub-NUT center, in which, unlike the above, both the 
metric and the dilaton go to zero as one approaches the singularity surface, and inside it the 
metric becomes negative and the space ceases to be Riemannian. The former is a common 
difficulty, while the latter singularity is the one inherited from the "one-loop approximation" 
of the Atiyah-Hitchin metric without taking into account the instanton effects. Nevertheless, 
despite the existence of these regions, we would like to emphasize that our construction is 
interesting in the following points: 

First, it clarifies the origin of orientifold-like objects in heterotic string theory; such objects 
are necessary to be included to neutralize the total NUT charge before T-duality. Indeed, 
their T-duals are very similar to orientifolds; the Z2 projection naturally arises from the bolt 
of the Atiyah-Hitchin space, and the ADM energy per unit area is negative. The second point 
we wish to make is that our 5-brane system is not just being a singular geometry, but suggests 
how the nonperturbative effects (by vortices in the corresponding gauge theory) correct the 
geometry. Although incomplete, we will consider in a later section the relevant (4, 4) sigma 
models whose leading-order moduli spaces are singular 5-brane geometries above but vortex 
corrections may change them. The problem amounts to finding a possible deformation 
of hyper-Kahler geometry with torsion, and the sigma model [13] analysis may provide a 
systematic way. 

The infinite periodic array of Taub-NUT centers that we describe is similar to three- dimen- 
sional ionic crystals, such as sodium chloride (NaCl), cesium chloride (CsCl), etc., where 
the Atiyah-Hitchin's and Taub-NUT's correspond to minus and plus ions in the crystal, and 
the Gibbons-Hawking potential plays the role of the Coulomb static potential of the ions. 
Generically, the Gibbons-Hawking potential for such a three-dimensional periodic array does 
not converge, even if the total NUT charge of a unit cell is neutralized. As we will see shortly, 
the condition for convergence is that the positive- and negative-charge centers are arranged 
in a unit cell in such a way that the leading difference operator of the multipole expansion 
of the potential becomes proportional to the discrete Laplace operator. This condition con- 



strains the locations of the Taub-NUT centers relative to the Atiyah-Hitchin centers, that is, 
the ones with negative NUT charge. Another constraint for the positions of the Taub-NUT 
centers is that they must be inversion symmetric with respect to every Atiyah-Hitchin bolt 
locus so that they are consistent with the Z2 identification around the bolt. Thus, while 
a real crystal is mathematically described as a three-dimensional torus, the "Taub-NUT 
crystal" is a certain orbifold. 

The plan of this paper is as follows: In section 2 we review the basic facts about the Atiyah- 
Hitchin metric. In section 3 we give a definition of our set up by introducing a three- 
dimensional periodic array of Taub-NUT centers, and examine the convergence condition 
for the Gibbons-Hawking potential. A relation to the K^ orientifold is discussed, and the 
classification of the crystal structure is also described there. In section 4 we take T-duality 
of the Taub-NUT crystal to obtain a 5-brane system with compact transverse dimensions, 
and then we embed it into heterotic string theory. At the end of this section we briefly 
discuss possible corrections to the singularity from the point of view of the nonlinear sigma 
model approach. Section 5 is devoted to a summary. An appendix contains some details 
of the Darboux-Halphen system, where comparisons are also made between the solutions in 
the literature. 



2 The Atiyah-Hitchin metric 

Let us first recall some basic properties of the Atiyah-Hitchin metric. The Atiyah-Hitchin 
metric is a smooth metric on a Euclidean four-manifold with a(n) ( ant i-) self- dual Riemann 
curvature without triholomorphic U{1) isometry. It is written as a special Bianichi IX metric 
of the form 

ds^ = a%^c^dt^ + a^al + feV^ + cV^, (2.1) 

where a, h and c are functions of t only, ctj's {i = 1,2,3) are the Maurer-Cartan 1-forms of 
SU{2) given by 

o"i = — sin ipdO + sin 6 cos iljd(j), 

o"2 = cos ip d6 + sin 6 sin iljd(l), (2.2) 

0-3 = dip + cos 9d(f). 

In terms of a new set of functions 

Wi = be, W2 = ca, W3 = ab, (2.3) 



a (sufficient) condition for a(n) (anti-)self curvature is that they fulfill the Darboux-Halphen 
system 

Wi+W2 = 2wiW2, 

W2+W3 = 2W2W3, (2.4) 

W3 + Wi = 2wsWi. 

For our purpose, a convenient representation of solutions is the one in terms of elliptic theta 
functions [Tit [15] 

Wj{t) = 2^1og^,+i(0,mt) (j = 1,2,3), (2.5) 

where a > is an arbitrary constant parameter. We will see that this a is related to the 
NUT charge of the asymptotic Taub-NUT metric. The obtained self-dual metric is thus 



ds^ = 2 JJ —log'dj{0,iat)]Udt'^ + 



dt " '^ ' ' )\ (|log^2(0,mt))2 

(2.6) 



Vi=3,4,2 / \ ^dt 

2 
(To Cr: 



(|log^3(0,2at))2 (|logt94(0,^at))^ 

In this representation of the solutions, Atiyah-Hitchin's original solution corresponds to 
a = TT, whereas Gibbons-Manton's solution to a = 27r (with replacing the radial coordinate 
t by —7]. See Appendix for more detail.). 

It is well known [16j (also see [17] for a review) that, on one hand, the Atiyah-Hitchin metric 
develops a removable "bolt" singularity at the origin (in a certain coordinate system), and 
on the other hand, it rapidly approaches to a Taub-NUT metric with negative NUT charge 
asymptotically. Although these are well-known facts, we include a brief derivation of them 
because they are essential in constructing a crystal structure in the space of Taub-NUT 
centers. These behaviors can be conveniently confirmed in our representation of solutions in 
terms of theta functions 



(i) The behavior near t = 00 .■ The bolt 

The elliptic theta functions are readily given as an infinite series as 

^3(0, r) = 5^g^"^ 

neZ 

MO,r) = ^(-1)M"', (2.7) 



neZ 



MO,r) = 5^g^("-^)', 



neZ 



with q = e^'^" . Putting r = iat, the t — )■ oo hmit of (12. 5 p can be obtained by taking the 
leading terms of the expansions: 

W2{t) *~°° 2^1og(l + 2e-™* + ---) 

= -47ra e"™* + ■ ■ ■ , 

w,{t) '^^ 2|log(l-2e~™* + ---) 

= +47ra e"™* + ■ ■ • , (2.8) 

«;i(t) *^~ 2|log(2e-f"* + ---) 

TT 
= --«+■■■. 

Therefore, a,b and c are 

a = ~ oZnae + ■ • • , 

2 W3W1 t-^00 vr 
b = ~ -aH , 

W2 2 

2 W1W2 t->oo TT 

c = ~ —a + ■ ■ • . 

W3 2 

The metric takes the form 

c/s^ *^°° 87ra(7rVe-2™*rft2 + 4e-2™V2) + ^a(a2 + a|). (2.9) 

This is a bolt; by a change of coordinates m = e^'^"*, the metric becomes 

ds'^ = 87Ta{du^ + Au'^cxl) + ^a{al + al). (2.10) 

Defining the new angle coordinates [16] such that 



0"! = d'?/' + COS ( 

(72 = — sin ■?/'(i6' + sin ^ COS 'j/'fi^, (2-11) 

(73 = cos ^d6' + sin 6 sin 7/'c?0, 

we see that this metric is regular at -u = (t = 00) ii < ip < it, that is, the points at the 
angles ip and ip + n must be identified. This Z2 transformation acts on the original angle 
variables as 

e ^ TT-e, 

H^ + 7r, (2.12) 

ijj I—)- —ijj. 



This means that the space around the Atiyah-Hitchin bolt locus gets Z2 orbifolded. Later 
this fact places a significant constraint on our construction of the Taub-NUT crystal. 

(a) The behavior near t = 0: The negative Taub-NUT 

To investigate the asympotics near t = 0, one needs to consider the modular S transforma- 
tions of the theta functions [T5l: 
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y/ at 1^3(0, iat) 
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tat J 


y/at 1)2(0, iat) 
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tat J 


Vat ??4(0,iat) 


Using (12.13p. we can similarly obtain 




W2{t) 


- 4- 


/ 1 i 
V Vat at 


Ws{t) 


= 4'°s 

dt 


( 1 i 
VVat at 




2at2 
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Wi{t) 
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/ 1 i 
Wat at 

•> 



(2.13) 



(2.14) 



where the ellipses stand for the terms containing a power of e <»* . If those contributions are 
discarded, the a, h and c functions reduce to 

2 t^o TT — 2at 



2at2 ' 
,2 t^o TT — 2Q;t 



'^ ^ ^sF- <^^^^) 



2 t-^0 

c ~ 



2a 



TT — 2at 
The metric approaches to 

2at4 7r-2at ^ 2at2 ^ ^ 2; 



-1 



2a \ / /2a 

1-^ (rfr^ + rV? + a2^))+f 1-^1 a3^ (2.16) 



where we have changed variables from t to r = {\/—t) ^ in the second hne. (I2.16P is a 



2a 



Taub-NUT metric with negative NUT charge —\-^ [16]. 

As was noted in p^, this approximate metric (I2.16p is singular at r = y^; this is an artifact 
of the absence of the exponential terms ignored in the theta functions. On the other hand, 
the smooth, self-dual metric ( 12. 6 p rapidly approaches to the approximate metric (I2.16P as 
r — )■ oo, and for an observer outside a radius of the NUT charge scale, it looks as if there is a 
Taub-NUT with negative NUT charge, which has an f/(l) isometry. Later we take T-duality 
of this approximate metric along this isometry. 



3 Taub-NUT crystal 

3.1 Infinite periodic array of Taub-NUT centers 

In the previous section we have seen that the Atiyah-Hitchin space can be approximately 
identified as a negative- charge Taub-NUT space at large distances compared to the NUT 
charge scale. Motivated by this observation, in this section we consider the Gibbons-Hawking 
metric for an infinite periodic array of Taub-NUT centers, where some of them in the unit 
cell have negative NUT charge. The idea is that, by periodically identifying such an infinite 
array, one obtains an S^ fibered three-dimensional torus with (anti-)self-dual metric, and by 
T-duality it is mapped to an NS5-brane system with compact transverse dimensions, where 
some of the 5-branes have negative tension. 

The Gibbons-Hawking metric [8] 

ds^Q^ = V{x)dx'^ + V^^{x)(dxQ + u ■ dx^, (3.1) 

V{x) = e + J^j^^^, 3.2 

VV = V X w, (3.3) 

is an anti-self-dual metrical in a four-dimensional Euclidean space with coordinates (x, Xg), 
X = {x,y,z). We have used the subscript 9, anticipating that this metric will be embedded 



^ Whether a metric is self-dual or anti-self-dual is a matter of convention. For instance, it changes if the 
ordering of the coordinates is changed from (x, xg) to (xg, x). We have chosen the sign of uj in such a way 
that xg is naturally identified as ip (up to a scale, without a sign flip) of the asymptotic metric (|2.16p . Also, 

the change of coordinates 1 1— > r = {J^t)~^ in (I2.16P makes a self-dual metric into an anti-self-dual one. 
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in string theory, e and m are constants. If e 7^ 0, the metric describes a multi-center Taub- 
NUT space, while if e = 0, the metric is the one for the A^^i singularity [18]. If £ 7^ 0, it 
can be arranged to 1 by making a suitable change of coordinates, so we can set e = 1 for the 
Taub-NUT case without loss of generality. Then if A^ = 1 and the center xi is taken at the 
origin 0, the potential V takes the form 

2m 



V{x) = l + {r= ^x'^ + y'^ + z^). (3.4) 



r 



Choosing the gauge 1^2 = 0, ci; can be determined (for the anti-self-dual case) as 

/ 2myz 2mxz \ 



I'hen the metric (13. ip reads 








rf4^(iv=i) = 


/ 2m\ 

1 + 

V r J 


1 dx' 


^ 1 + 2™ V 


By a change of coordinates 










X 


= 


r sin 9 cos 0, 




y 


= 


r sin 9 sin 0, 




z 


= 


r cos 6*, 




Xg 


= 


2mV', 



(3.6) 



(3.7) 



(13. 6p is reduced to a Taub-NUT metric with NUT charge 2m: 

dslH{N = l) = fl + ^ydr' + r\al + al)) +4m' (l + ^] \l (3.^ 



Going back to the general N case, we impose the following conditions on the configuration 
of centers: 

1. For every location of a negative NUT charge Xq, the configuration is invariant under the 
Z2 transformation : 

X — xq I— !■ — (x — xo) (3.9) 

(called inversion in the terminology of crystallography). 

2. The superposition of infinitely many Gibbons-Hawking potentials converges. 

The first condition is imposed in order for each negative- charge center to consistently repre- 
sent the asymptotic region of the Atiyah-Hitchin geometry. To recover correctly the original 



Atiyah-Hitchin metric in the non-compact limit where the negative-charge NUT is isolated 
from the positive ones, the image of the inversion (13. 9 p needs to be identified, as we saw in 
the previous section. 

The second condition is obvious; a Gibbons-Hawking potential from a single nut falls like 
l/|x|, and due to the condition 1 the leading asymptotic behavior of the superposition 
is l/|a;|'^ generically, and therefore the summation over a three-dimensional lattice diverges 
logarithmically. The superposition of the potentials converges if and only if the positive- and 
negative-charge centers are arranged in a unit cell in such a way that the leading order term 
of the multipole expansion vanishes. We will present a precise description of the condition 
in the next section. Note that the situation is completely analogous to the Coulomb static 
potential in a crystal system, and any actual ionic crystals consisting of two oppositely 
charged ions will satisfy this condition. 

In addition, we also assume for simplicity that 

3. There are no orbifold fixed points besides the locations of negative charges. 

Then a little thought shows that the Gibbons-Hawking potential V{x) takes the form 

v{x) = 1+ Yl E ^^"'"'^'Ha^-E-=i^,y,), (3.10) 

ni,n2,n3£Z <:i,i:2,e36{0, 2 } 






2m 



^ - E,=i ^jVj 



k=l 



X ~ X 



{<^l,<^2,>^3) 



V^3 



2m 



X + Xf^ — 2^j=i ^jVj 



(3.11) 



where yj G R^ (j = 1, 2, 3) are the period vectors of the array. 



The unit cell of the periodic array is a parallelepiped specified by y^ {j = 1,2,3), which 
contains eight (e^ G {0, |} for j = 1,2,3) subunits (Figure). Each subunit specified by 
(^1; ^2, £3) of the array consists of one negative-charge center with NUT charge —4:mN^'^'^'^^''^^' 
and A^('^i''^2,e3) pairs of positive-charge centers with a common NUT charge 2m, which are 
placed symmetrically with respect to the position of the negative-charge center. The relative 



(<:i,<:2,<:3) 



positions of the fcth pair (A; = 1, . . . , , A^('=i''=2,<:3)^ ^re ±x 
arises from these centers belonging to the subunit (ei, €2, 63). 



The potential ^(^I'^a-^s), 



x) 



.^B^ 





Figure 1: (a) A Unit cell and (b) a subunit of the Taub-NUT crystal. The dots are the 
positions of negative-charge centers, of which the circled ones are identified periodically. 
Positive-charge centers, which are not depicted in these figures, must be placed symmetrically 
with respect to the inversion through any dot. 



We now periodically identify points as 






(3.12) 



and also take an orbifold 



[X, Xg) 



-x + 2 ^{rij + ej)yj, -Xg | {rij e Z, e^ e {0, i}), (3.13) 



as we mentioned above. The metric 03.11) with (I3.10p . (l3.1ip and (13. 3p is invariant under the 
operations (I3.12p and (I3.13P and the identifications can consistently be done. 



Note that if there were only a single negative-charge center, (I3.12p would identify points 
in spaces with opposite orientations and the quotient would become non-orientable. Also, 
there would appear orbifold fixed points besides their own locations. This is why we have 
introduced eight negative-charge centers in the unit cell. 



3.2 Convergence condition 



The metric (13. ip with (13.100 . (13. lip and (13.30 is an infinite, periodic generalization of the 
Gibbons-Hawking metric. In this section we consider in what circumstances the potential V 
converges. If the positive and negative centers are in generic positions, the potential (13. lip 
falls like l/|if|^ since it can be written as a sum of second order differences, and therefore the 
summation over a three-dimensional lattice diverges logarithmically. However, the exception 
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to this is where the leading order term of the multipole expansion of 

J2 V^'''''''-'\x) (3.14) 

vanishes, that is, the positive- and negative-charge centers are arranged in a unit cell in such 
a way that the leading difference operator of the multipole expansion is proportional to the 
discrete Laplace operator. 

The simplest example of a convergent potential is the case where the lattice is a cubic lattice 

y, = 2an, (j = 1,2,3), (3.15) 

ni = (1, 0, 0), n2 = (0, 1, 0), n3 = (0, 0, 1), (3.16) 

and all V^'^^''^^''^^\x)^s are obtained as translations of an identical potential: 

VReoM^) ^ 2m(-^ + ^/' ^ -i + i-^A -I IV (3.18) 

\ \x\ ^—^\\x — danj\ \x + oarij] J I 

where a is half the period of the array and < 5 < 1 is a parameter (Figure |2] (a-1), (a-2)). 
I 

Indeed, we can rewrite (I3.18P as 

2m| --^ + V ( ,^ \ ^ . + .^ ^^ 
\ \x\ ^-^ \\x — oarijl \x + dan 

3 

= 2m^ 



i=i ^'"" '""■" '" ■ '"'^'^^ 



, ^ ,x + SanA \x\ J \\x\ \x — SanA 



^ 1 

2m(5V^(rr,-V)2p- + 0( 



a^) 



i=i 



= 2m(5VA^ + 0(a^), (3.19) 

\x\ 

where A is the flat space Laplacian. The first term vanishes, and on dimensional grounds 
the rest falls like |x|~^ at least. Therefore the summation over the lattice converges. 



^The actual crystal of rhenium trioxide (ReOs) realizes this configuration for S — 1. Strictly speaking, if 
the Taub-NUT positions change from S y^ 1 to 6 — I and two Taub-NUT centers come on top of each other, 
then the periodicity of the xg direction must be changed so as to avoid the NUT singularity. Here, we are 
interested in the convergence of the three-dimensional potential and ignore this complication. 
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(a-1) 



(a-2) 




(b-1) 



(b-2) 



Figure 2: Examples of convergent configurations. The unit cells of (a-1) the ReO^-tjpe 
and (b-1) the CsCl-tjpe crystals are depicted, and (a-2), (b-2) are their sub-units, respec- 
tively. The blue spheres show the positions of the negative-charge ( "Atiyah-Hitchin" ) centers, 
whereas the red ones are the positive-charge ( "Taub-NUT" ) centers. 



Another simple example of a convergent potential is the case in which the lattice is again 
cubic and a subunit has four (pairs of) Taub-NUT's per single Atiyah-Hitchin (Figure [2] 
(b-1), (b-2)): 



Vcsciix,6) 



Vcsciix - Y.^j^l'^jVj, S), 



2m 



\x\ ^-^ 



a,6,€3=±i F - '^« Yl,i=\ Ci^j 



(3.20) 
(3.21) 



More generally, if one requires that for each subunit (ei,e2,e3) the leading-order term of 
^(ei,e2,e3)('^^ lu the large- |x| expansion vanishes, the condition for the relative positions of 
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the Taub-NUT centers is that the 3 x A^'^'^^-'^^'^^^-niatrix 

I^(ei,e2,e3) — f d^l'^2,e3) -j(ei,e2,<:3) -.(ei.ea.ea) A (o 99^, 

should satisfy 

js^(.i,.2,.3)js^(^i,^2,e3)^ = ^jg (323) 

for some constant /c, where Is is the 3x3 identity matrix. It is easy to see that the two 
examples above fulfill this condition. 

Note that (I3.23P is a sufficient condition for V{x) to converge; it is enough if the sum (13.141) 
over the unit has the asymptotic fall-off is not slower than |a;|~^, and each V^'^^''^'^''^^\x) need 
not subunit-wise. Of course, in such a general case the crystal has less discrete symmetries, 
however. 

3.3 Taub-NUT crystal as an approximation of K2, 

The negative-charge metric (12.161) is singular, but adding exponential correction terms the 
singularity disappears and the metric becomes completely smooth to be the Atiyah-Hitchin 
metric. These exponential corrections are regarded as coming from certain instantons |T9l 
[20l [TSj . In the language of the corresponding gauge theory, this is a strong-coupling effect 
of the non-abelian (= SU{2)) gauge theory at low energies, where the location of the singu- 
larity surface corresponds to Aqcd- It is reasonable to expect that there will also be such 
corrections even after the three-dimensional lattice is compactified to be a torus. If this is 
true, and if these corrections can smooth out the singularity as they do in the non-compact 
case, then the resulting space must be a K3 surface, as it is the only nontrivial smooth 
compact hyper-Kahler four-manifold. The compactified periodic array can then be regarded 
as an approximation of a K^ surface (which is only valid except near the locations of the 
negative-charge centers). However, this can only happen if 

<:i,e2,e3G{0,|} 

in a unit cell since the Euler number of K^ is 24. This means that (apparently, including 
their Z2 images) there are 32 positive-charge Taub-NUT centers in a unit cell before the Z2 
identification. Note that although the number of negative-charge centers is apparently eight, 
it must be halved when counting the number of bolts due to the Z2 identification. (This is 
the same as the Euler- number counting of the K^ orbifold [21j , where the number of excised 
points is not equal to the number of the fixed points itself, but it must be divided by the 
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order of the orbifold.) Therefore, since a single nut adds one to the Euler number while a 
bolt adds two to it, the total Euler number amounts to 

1 X 16 + 2 X - = 24. (3.25) 

In this case, the compactified and orbifolded periodic array is related by a chain of dualities 
to the K^ orientifold, where the sixteen Taub-NUT's here correspond to the same number 
of D9-branes. 



3.4 Classification of the crystal structure 

Crystals are known to be classified by their symmetries. The crystal structure is characterized 
by the following two aspects: One is the lattice that determines its translational symmetries, 
and the other is the fundamental structure whose periodical repetitions generated by the 
lattice yields the whole crystal (Figure |3]). Three-dimensional crystals fall into seven classes, 
the crystal systems, depending the shape of the parallelepiped of the unit cell that the crystal 
lattice specifies. Four of the seven crystal systems, the monoclinic, orthorhombic, tetragonal 
and cubic systems, are further divided into several sub-classes, as shown in Table [TJ The 
definitions of the edge lengths and angles of the parallelepiped are shown in Figure HI 



Two crystal lattices belonging to the same crystal system but different sub-classes (that is, 
two ones in the same row in Table [T]) possess the same point group symmetries (that is, 
the group of symmetry operations that leave at least one lattice point fixed) but different 
translational symmetries. These, in all, fourteen crystal lattices are known as the Bravais 
lattices in crystallography. The classification into the fourteen Bravais lattices is further 
refined into 230 crystal types according to their discrete symmetries called space groups, 
which are generated by the point group and the group of translational symmetries. 




+ o/ = 




1 attire fundamental crystal 

Id L Lice structure structure 



Figure 3: The crystal structure. 
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crystal system 



{length, angles} 
{symmetry element} 



Bravais lattice 



C 



triclinic 



{a 7^ 6 7^ c, a 7^ /3 7^ 7} 
{no symmetry element} 



mono clinic 



a 7^ 6 7^ c 

a = 7 = 7r/2,/3> 7r/2 

{ 1-diad} 





orthorhombic 



a 7^ 6 7^ c 

q; = /3 = 7 = 7r/2 

{mutually orthogonal 3-diad} 





hexagonal 



a = h ^ c 
7 = 27r/3 
{1-hexad} 




trigonal 



a = h = c 

7 7^ 7r/3,7r/2,cos"^(-|^ 
{1-triad} 



tetragonal 



a = h ^ c 

a = f3 = J = 7r/2 
{1-tetrad} 




cubic 



a = b = c 

a = (3 = J = 7r/2 
{4-triad} 





Table 1: The Bravais lattices. 
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Figure 4: The definitions of tlie edge lengtlis and angles in Table [H 



One of the important aspects of such discrete symmetries of a crystal is that they might be 
used to explain an origin of discrete flavor symmetries |22] in nature (if any), though in the 
present case in the compactification to six dimensions. Indeed, it has been known for a long 
time that there are localized p3j 30 hypermultiplets on each symmetric 5-brane (obtained by 
T-duality in the next section) , whose existence is ensured by anomaly cancellation [211 125] , 
and these hypermultiplets are to be observed as flavors in the warped compactification. It 
would be interesting to generalize the argument to a similar periodic array of intersecting 
5-branes. 



Taub-NUT crystals with a highly symmetric structure are very limited. For instance, let us 
consider cubic lattices with "triads" (that is, symmetry axes of rotations by 120 degrees) 
that have 16 positive-charge Taub-NUT's (Euler number = 24) as in section 13.31 Then it 
turns out that the possible type is either of the three shown in Figure [51 




Figure 5: Cubic lattices with triads that have 16 positive-charge Taub-NUT's. (a) and (b) 
are of the P type and their subunits are shown there, (c) is of the I type and hence a region 
larger than a sub unit is depicted so that its symmetry may be manifest. 
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4 T-duality to NS5-brane system 

4.1 0(1, 1) dual and embedding to heterotic string theory 

It is also a well-known fact that the NS5-brane [23l [26] and the Taub-NUT space are related 
by T-duality [271 [2H], which we briefly review here. The NS-NS sector Lagrangian is given 
in the string frame as 

2k'C^sns = Ee-"'(R + AdM^d''<^-^HMNpH''''''Y (4.1) 

where E is the determinant of the vielbein Ej^, related with the string- frame metric Gmn 
(M, iV = 0, . . . , 9) as Gmn = Ej^fE^rjAB- Using the Lorentz invariance, one can cast the 
vielbein into the form 

( E- A,E\ 
Ej^ = I y (4.2) 

where /z (a) is a nine-dimensional curved (flat) index. If the metric and all other flelds do 
not depend on x^, then decomposing the vielbein as the Lagrangian can be reduced to an 
0(1, 1) sigma model coupled to gravity, an B fleld, two abelian vector flelds and a dilaton, 
which is invariant under the following transformation [29j : 



(4.3) 



771 new a 


= 


E,-, 


/I new 


= 


Bfj.9, 


pnew 


= 


E-\ 


nnew 


= 


B^u - 2 A [,^5^] 9 


onew 


= 


Afj.} 


^new 


= 


$ - log E. 



Therefore, starting from a solution of the equations of motion, one can obtain another new 
solution by applying this transformation to the old one. 

We apply the T-duality transformation to the Gibbons-Hawking multi-center metric [S]. If 
the Gibbons-Hawking metric (13. ip is embedded in eleven dimensions and the Xg direction is 
compactifled on S^, the metric describes A^ D6-branes in type HA theory; this is a string 
theory analogue of the GPS monopole. Alternatively, we embed (13. ip as a four-dimensional 
metric in ten dimensions so that the vielbein becomes 

le 

E^ = \ Vhs V'"^Q I , (4.4) 

\/-5 
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while other NS-NS fields are set to zero. Then the 9 + 1 decomposition reads 



V 





3 

(/i = 0,...,5), 



vh 



\ 



n 



n 



^x (yU = 6), 

bJy (^ = 7), 

w^ (/i = 8), 
E = V-^ (4.5) 

B^,y = 0, 

B^9 = 0, 
$ = 0. 

Using the transformation rules (14. 3p . we obtain the T-dual configurations as 

ic n 

rpnew a I 

'3 



^ \ \/n. 



4 new _ 


0, 






rpnew 


v^^ 






onew 
-D fiu — 


0, 






rjnew 
-D m9 — 


' 

Wx 

< 




= 0, 
= 6) 

= 7) 
= 8) 


^new ^ 


^logK 







(4.6) 



5), 



or 



'^"^new = rfijdx^dx^ + V{x)S^^dx^dx'^, (4.7) 

^;:^^ = eV9^p^(^), (4.8) 

g2$- ^ v{x), (4.9) 

where e^'^^g is normalized to take values ±1. (14. 71) - 04. 9 p describe A^ smeared NS5-branes. 

It is now obvious that, to obtain the T-dual of the infinite periodic array that we described 
in section 3, we have only to replace V with (IS.lOp . (13. lip in (I4.7p - (l4.9p . as the T-dual 
transformation does not depend on the explicit form of V{x). Thus we obtain a three- 
dimensional periodic array of, in this case, parallel NS5-branes with one (x^) dimension 
smeared and compactified on a circle. The Z2 projection (13.130 must also be imposed. 



It can also be easily verified that ^^^0 ^ HJ^a^, where fij^^^ is the spin connection derived 
from E^'^'"^, belong to different SU{2) subalgebras of the transverse space holonomy group 
5*0(4). Therefore, the 5-brane system has a hyper-Kahler with torsion geometry, and hence 
preserves one half the supersymmetries. We can then construct in the standard fashion a 
supersymmetric heterotic background by identifying some part of the gauge connection to 
be equal to the generalized spin connection: 

fi°r^±if-^ = (A^U, (4.10) 

where a, (3 on the right hand side are the gauge indices specifying some SU{2) subalgebra 
of Eg X Eg or 50(32). 



4.2 T-dual of negative- charge Taub-NUT as heterotic orientifold 

The Z2 identification through every location of the "Atiyah-Hitchin" brane shows that they 
are orientifold-like objects in heterotic string theory. Another property which supports the 
identification is their tensions: In the string frame, they have negative ADM energy per unit 
area as can be easily verified. Of course, this is made possible because they have singularities. 

In order to construct a 5-brane system with compact transverse dimensions, we first con- 
sidered an infinite periodic array of Taub-NUT centers. To make the Gibbons-Hawking 
potential converge, we were naturally led to introduce negative charge centers as a part of 
it. Then by T-duality, they turned into branes with negative tension. That is, they need 
to appear, in this construction, for the transverse space to be compactified. This may be 
viewed as in accordance with the statement of [1]. 



Also, in [251 EO], a brane compactification of the E^ x E^ heterotic string theory was con- 
sidered using smeared intersecting 5-branes and the spectrum of fermionic zeromodes was 
computed. There, some of the intersecting branes were required to have negative tension for 
a convergent harmonic function of the metric. However, without a microscopic definition of 
orientifolds in heterotic string, their origin was obscure. The T-dual of (the asymptotic form 
of) the Atiyah-Hitchin space thus may provide an understanding of such negative tension 
branes in heterotic string theory (despite the singularity; see next section). 
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4.3 Discussion: Vortex corrections to the singularity 

So far we have considered a system of brane compactification of string theory in the frame- 
work of supergravity. In any case, it is not useful when the string couphng is not small, 
or when the fields vary rapidly compared to the string scale. In this section we discuss the 
possibility of probing the geometry near the 5-branes by using the sigma-model approach. 
Of course, heterotic string has no D-brane probes, but one may use them first in type II 
theories as a systematic means of obtaining a hyper-Kahler geometry (with torsion), and 
then it may be converted into a heterotic background by the standard embedding. 

The three-dimensional A^ = 4, SU{2) supersymmetric gauge theory obtained by a dimen- 
sional reduction of the four-dimensional N = 2 gauge theory without hypermultiplets has the 
Atiyah-Hitchin space as its moduli space [191 [20]. In three dimensions, the U{1) gauge field 
is dualized to a scalar, so that there are four real scalars whose expectation values parame- 
terize the Coulomb branch moduli space, which must be some hyper-Kahler manifold [31]. 
To one-loop, the metric is nothing but the negative-charge Taub-NUT that we have used in 
the Gibbons-Hawking metric. By taking into account the corrections due to monopoles, it 
was shown that the smooth Atiyah-Hitchin space was singled out from other hyper-Kahler 
manifolds [121 [20] • Originally, the Atiyah-Hitchin space was discovered [D] as the moduli 
space of slowly-moving two monopoles in a four- dimensional gauge theory [32]. The appear- 
ance of the same moduli space of different gauge theories was explained [33] by considering a 
D3-NS5-brane system, where the two gauge theories arise as low-energy theories on different 
branes. 

Now suppose that we consider a further dimensional reduction to two dimensions. This is a 
two-dimensional (4, 4) simga model, consisting of a hypermultiplet and a twisted hypermul- 
tiplet, both in the adjoint of SU{2). In two dimensions the vector field have no propagating 
degrees of freedom, and there is no need to dualize anything. The degree of freedom that 
the three-dimensional dual scalar had is now carried by one of the scalar components of the 
twisted hypermultiplet. This is essentially the Buscher's T-duality, and the sigma model is 
thus expected to describe the T-dual geometry of (in the SU{2) case) the Atiyah-Hitchin 
space. 



This set up was already considered long time ago by Diaconescu and Seiberg [51]. Indeed, 
in the SU{2) case, they obtained a one-loop (exact) sigma model metric precisely the same 
as the singular negative-tension brane that we obtained as the T-dual of the negative-charge 
Taub-NUT. (See also [35].) The relative minus sign in the harmonic function is a reflection of 
the non-abelian nature of the gauge theory, and in the f/(l) case the relative sign is positive. 
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and that corresponds to the ordinary Taub-NUT space. However, unhke in three dimensions 
where monopoles can exist, there are, at least in the naive sense, no further non-perturbative 
contribution from vortices in two dimensions on this Coulomb branch; there is no vortices 
in the reduced two-dimensional theory of the pure SU{2) gauge theory broken to U{1) since 
7[i{SU{2)/U{l)) is trivial. Still, we note that it has also been known for some time that, in 
a sense described in |H6] , the vortex moduli space of the Hitchin system is well-defined and 
its volume can be computed 



Realizing the hyper-Kahler geometry with torsion on the Higgs branch is also interesting. 
In |37], the smeared 5-brane metric was realized as the Higgs branch moduli space of a U{1) 
(4, 4) sigma model, and it was argued that the world sheet instantons (= vortices) correct 
the metric into the Poisson re-summed form of the localized 5-brane metric compactified on 
a circle. It would be interesting if one could generalize this to the negative-tension/Atiyah- 
Hitchin case. The hyper-Kahler quotient construction for the singular two-monopole moduli 
space (= the Atiyah-Hitchin space without instanton corrections) is already known [38J. 
What puzzles us about the result is that, in that case, one needs to start from a flat hyper- 
Kahler space with an indefinite metric, which would mean that some of the scalars of the 
sigma model would have to have a wrong-sign kinetic term. The realization on the Higgs 
branch requires further study. 



5 Summary 

In this paper we have taken a first step toward realizing warped compactifications using 
NS5-branes in heterotic string theory. We have first considered the Gibbons-Hawking metric 
for an infinite periodic array of Taub-NUT's, and have taken T-duality to obtain a parallel 5- 
brane system with a compact transverse space having the SU{2) structure, which allowed us 
to obtain a supersymmetric heterotic 5-brane system via the standard embedding. In order 
for the Gibbons- Hawking potential to converge, there have been included Taub-NUT centers 
having a negative NUT charge, which have been identified as the asymptotic forms of the 
Atiyah-Hitchin metric. This identification has forced us to make a Z2 projection through the 
positions of the negative Taub-NUT centers. We have described the convergence condition 
for the locations of the centers. We have shown that in order for the periodic array to be an 
approximation of K3, the number of the positive charge Taub-NUT's must be 16 in a unit 
cell. We have seen some parallels between the periodic array of the Taub-NUT centers and 
ionic crystals, and structures of the "Taub-NUT crystal" can be similarly classified in terms 
of the Bravais lattices, of which we have given a brief review. 
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Since we have used the asymptotic form of the Atiyah-Hitchin metric as a part of the con- 
struction, the metric (13. ip with fl3.10p .f l3TTT]) and fIS.Sp have singularities near the locations of 
the negative charges [HI [12], and it only makes sense at distances larger than the charge scale 
from those charges. Despite the singularities, the compactified/orbifolded infinite periodic 
array has significance in that (i) it clarifies the origin of negative tension branes in heterotic 
string theory, and (ii) it basically gives the leading-order behavior of such objects whose cor- 
rections from nonperturbative effects might be systematically analyzed by the sigma model 
approach as we discussed in the text. 



Appendix The Darboux-Halphen system 

A.l The Bianchi IX self-dual metrics 

The Bianichi IX metric: 

ds"^ = a%^c^dt^ + a^al + 6V^ + cV^. (A.l) 

a, 6 and c are functions of t only, ctj's {i = 1, 2, 3) are the Maurer-Cartan 1-forms of SU{2) 
given by 

o"i = — sin ipd6 + sin 6 cos iljd(j), 

a2 = costpdO + sin 6 sin ipdcj), (A. 2) 

(73 = dip + cos ( 



If one requires that the spin connection co^ap is self-dual 



then a, 6 and c must satisfy 



a 


1 ^ -rS 

a'-b'- c2 

1 


abc 
6 


' 26c ' 
62 - c2 - a2 


abc 
c 


' 2ca 

c^-a'- 62 

1 


abc 


' 2a6 ' 



(A.3) 



(A.4) 



where ' = ^- These conditions are satisfied by the Eguchi-Hanson metric written in the 
form (lA.ip . More generally, if one requires that the Riemann curvature tensor is self-dual: 
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the conditions are [T6 



a a? — l)^ — c^ 
abc 2bc 

b b"^ — c^ — a"^ 

abc 2ca 

c c^ — a^ — b"^ 

abc 2ab 

with 

Since (lA.T^ is invariant under the sign flips 



a 


t-)- 


+a, 


b 


I-)- 


-b, 


c 


I-)- 


-c, 


Xa 


I-)- 


+ Xa, 


Xb 


H-). 


—Xb, 


Xr 


I-)- 


-Xr 



(A.8) 

and analogous cyclic generalizations, the consistent values of Aj's are essentially (Aq, A^, Ac) = 
(0,0,0) or (1, 1, 1) only. The latter choice corresponds to the cases of the Taub-NUT and 
the Atiyah-Hitchin spaces. 

In the original Atiyah-Hitchin paper [9], the overall signs of all the three equations are 
flipped, which correspond to an aniz-self-dual curvature in our convention. See below for 
more detail on the Atiyah-Hitchin metric in the original literature. 



A. 2 The Darboux-Halphen system 

This review is based on [14J . As was shown in [9] , by a change of variables 

wi = be, W2 = ca, wz = ab, (A. 9) 

the equations flA.6p with A^ = Aj, = Ac = 1 reduce to 

Wi+W2 = 2wiW2, 

W2+W3 = 2w2Ws, (A. 10) 

ws + wi = 2wsWi, 
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which is known as the Darboux-Halphen system. As we mentioned in section 2, 

w,{t) = 2^1og^,+i(0,zat) (j = 1,2,3), (A.ll) 

{a > 0) is a solution p3| [T5] . This can be shown by using the composition formula of the 
theta functions and the fact that the theta functions fulfill the heat equation. 

The original Atiyah-Hitchin solution can be written in this form as 

wf^'iv) = -2^iog^^.^i(0,-z7rr/) (j = 1,2,3), (A.12) 

while the solution obtained by Gibbons-Manton [16] is 

wf^if]) = -2|-log^,+i(0,-2z7rr7) (j = 1,2,3). (A.13) 

As we remarked at the end of the last subsection, (]A.12p has extra minus signs since they 
satisfy the anti-self-dual equations (see flA.14p ). Comparisons of the solutions are made in 
the next section. 



A. 3 Comparison with solutions in the hterature 

In this section we derive eqs.( ]A.12|) and (lA.13p . 



A. 3.1 Atiyah-Hitchin's original solution 



Let us first consider (]A.12p . In [9], a set of solutions are given to the Darboux-Halphen 
system 

w[+W2 = —2WiW2, 

■^2 + "^3 = — 2U72W3, (A. 14) 

which corresponds to a, 6 and c satisfying 

a' a^ — b"^ — (? 

abc 2bc 

V_ _ b^-c^- a^ 

abc 2ca 

abc 2ab 
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where the prime ' denotes j-. The metric is assumed to have the form (jA.ip with the 'radial' 
coordinate rj replaced instead of t. The relations between w/s and a,b,c are flA.Qp . The 
solutions are given in terms of the coordinate /3 such that 

(A.16) 









rf/3 _ ^2 
drj 


where u = 


= u{P) is a 


function of /3 


satisfying 

d'^u u^ 
d(3^ ' 4sin2/3 ~ 


A solution to flA.17) 


is explicitly 





0. (A.17) 



m(/3) = (2 sin (3)^ K {sin ^), (A.18) 



where 

Jo (1 — K^^sm 0j2 
is the complete elliptic integral of the first kind. According to [9j, we set 

4H- du u^ ^ 

'"^" = -"Ji^ + f-^;. (A-20) 

(ip 2 smp 

ATT du v^ 1 

Wi = -u 



d(3 2 sin/3 



Using ( JA.181) . we have 



w^"" = — ^logir(sin^' 
^ dr] ^ ^ 2' 



AJ^ c? ,_ ( ^ (3 ^^,_^ /3, 



"^^ = -^i°s(^^^^|^^^^^t))' ^^-^i^ 



AH ^ i„^ / „„„/5 r^/.:„ Z^- 



The elliptic theta functions and the elliptic integrals are known to be related as 



TT 



Vr^K{k) = -r?4(0,r(A;))2, (A.22) 

kK{k) = 1^2(0, r(fc))^ 
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where the modulus T{k) is given by 



Let us consider W2^ . Using these equations, we can further rewrite the expression of w^^ 



as 



w^ = -2^log^3|0,.§^|. (A.24) 



di] \ Kisin 



2- 



Using the differential relations among the elliptic functions one can show the relation 

4r(sinf) = ^^ ^. (A.25) 

d(3 ^ ^' 2sin/3ir(sinf)2 ^ ' 

The right hand side is equal to -:^^- Comparing flA.25p with (IA.16p . we find 

^^(sin f ) . dr] 

If we choose the arbitrary integration constant to be 0, we have 

r(sinf) = -Tiir], (A.27) 

and 

w^" = -2—\og'd3(0,-i7iri). (A.28) 

dr] 

Expressions for w^^ and w^^ can be obtained similarly. 

A remark is in order here. Since the imaginary part of the elliptic modulus of must be 
positive, flA.27p implies that rj takes its value from — oo to 0, which is puzzling if it is 
regarded as a radial coordinate. Since the system of equations (lA.14p is invariant under the 
simultaneous sign flips Wj t— )■ —Wj and r] i— )■ — r/, one might wonder if 

wf = -2—\og^j+iiO,+i7rri) (A.29) 

also work. They do not, however, because, although (IA.29P still solves flA.Mp . a^, P and c^ 
derived from (1A.29P with (lA.Op turn out to be negative. This problem can be cured by either 
changing the relations ( JA.9P to ones with minus sign: wi = —be, etc., or replacing ( JA.14P 



with the self-dual equations ( JA.IOP without minus signs on the right hand sides; in any case 
one is forced to modify some of the assumptions made in |9]. We have numerically checked 
that (1A.12P (and not (JA.29P ) is correct by Mathematica. 



While, of course, these complications are not essential, we adopted f lA.lOp instead of f lA.14p 



in the text for a simpler notation with less minus signs. Note that, unlike rj, the coordinate 
t in (lA.lip takes its value from to +cx). 
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A. 3. 2 Gibbons-Manton's solution 



In [16], the radial coordinate was chosen so that the Bianchi IX metric took the following 

form: 

52 
ds'^ = —dr^ + a'^al + b'^al + c^al (A.30) 

a, b and c which solve 



r da 


1 r 2 2 

a —b — c 


b dr 


2bc 


r db 
b dr 


b'-c'-a^ ^ 
2ca 


r dc 
b dr 


c^-a'-b^ 
2ab 



(A.31) 



give an anti-self-dual metric. The solution was similarly given in terms of Wj = wf^^ (j 
1,2,3) flA9l) as 



w^^ = -fsin/3 



dr 

dp' 

dr f^ 



with 



wf^' = -fsin/3^ + -(l-cos/3), (A.32) 

w^^^ = -fsin/3— - — (l-cos/3) 

(Jjj:) Zi 



2is: ( sin I ) . (A.33) 



To derive ( lA.13p . we change variables so that the metric takes the form ( JA.ip : 



abc drj = — dr. (A. 34) 

f 



rj is the new radial coordinate. Then (IA.33P implies that 

dfj 1 



dl3 2u(/3)2' 

where u{f3) is the function (lA.lSp . On the other hand, 102^^^ can be similarly written as 



(A.35) 



,GM " 1_ _, r^ / _•__ P 



w^^ = -— logiTlsin^) (A.36) 



d , / ir(cosf) 

■ -log^3 0, V • J 
dr] \ K{sm f ) 



-2-^log^3|0,^ J . l: ]. (A.37) 
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Comparing (1A.25P with (IA.35p . we have in this case 



and therefore 



r ( sin - ) = -2mfi, (A.38) 



nu ^ _2^\og^J0,-2i7rfi). (A.39) 

dr] 



The proofs for w^ and Wi are also straightforward. 



We have shown that the value of a flA.lip for Gibbons-Manton's solution is 2tt. Plugging 



a = 2tt into fl2.16p . we find agreement with the asymptotic metric obtained in 
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